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FLATNESS FOR A STRONGLY DEGENERATE 1-D 
PARABOLIC EQUATION 

IVAN MOYANO 


Abstract. We consider the degenerate equation 
dtf{t,x) - dx {x°‘dxf) {t,x) = 0, 

on the unit interval x £ (0,1), in the strongly degenerate case a G [1, 2) 
with adapted boundary conditions at a; = 0 and boundary control at 
a: = 1. We use the flatness approach to construct explicit controls 
in some Gevrey classes steering the solution from any initial datum 
/o G L^{0, 1) to zero in any time T > 0. 

Keywords— partial differential equations; degenerate parabolic equation; 
boundary control; null-controllability; motion planning; flatness. 


1. Introduction 


We consider the following control system 


( 1 . 1 ) 


dtf{t, x) - f{t, x) = 0 , 

f(t,x)l^=o = 0 , 
f(t, 1 ) = u(t), 
f(0,x) = fo(x), 


(t,x) € (0,T) X (0,1), 
t € (0, T), 
t G (0, T), 

X G (0,1), 


where the state is the solution f{t,x) and the control is the function u{t). 
The parameter a G [1,2) is fixed through the whole article. 

The aim of this work is to construct explicit controls u for the null- 
controllability of system (1.1) in finite time T > 0, using the flatness method. 


1.1. Main result. We will make use of the Gevrey class of functions. 

DEFINITION 1.1. Let s G M"*" and ti,t 2 G M with ti < t 2 - A function 
h G '^°°([ti, t 2 ]) is said to be Gevrey of order s if 


3M, R > 0 such that 


sup 

t\<r<t2 




< 


M{niy 


We then write h G ?f^([ti, t 2 ])- 


Before stating the main result, we have to recall the notion of weak solu¬ 
tions of the inhomogeneous system (1.1). 

DEFINITION 1.2 (Weak solutions). Let /o G L^(0,1), T > 0 and u G 
H^{0, T). A weak solution of system (1.1) is a function f G ^^([0, T]; L^(0,1)) 

1 






2 


I. MOYANO 


such that for every t' G [0, T] and for every 

(1.2) 'if G ^1([0, t'y,L\o, 1)) n ^0([0, t'yH\o, 1)) 

such that 


(1.3) 
one has 


{x°'dx)ip{t,x)\x=o = ifit,!) =0, Vt G [0,t'], 


f [ f{t,x){dt'tf + dx{x'^dx'if)){t,x)dtdx 
Jo Jo 

f f{tyx)'ip{t',x)dx- f fo{x)'if{0,x)dx + [ u{t)dx'ip{t,l)dt. 
Jo Jo Jo 


As we show in Section 2 (see Corollary 2.2), system (1.1) has a unique 
weak solution under suitable assumptions. Our main result is the following. 


THEOREM 1.3. Let /o G L^(0,1), T > 0, r G (0,T) and s G (1,2). 
Then, there exists a flat output y G ?f^([r, T]) such that the control 


(1.4) u{t) 


0 , 

^poo __ 

2^k=0 (2-a)2'=fe!n,ti(l+ff^)’ 


ift G [0,r], 
ift G (r,r], 


steers to zero at time T the weak solution of system (1.1). Furthermore, the 
control u belongs to ?f^([0,T]). 


1.2. Previous work. 


1.2.1. Null-controllability. The null-controllability of system 

dtf{t,x)-dx{x°‘dx)f{t,x) = lu,{x)v{t,x), (t,a;) G (0,r) X (0, 1), 

(x“5^)/(t,x)|a;=o = 0, tG(0,T), 

/(t,i) = o, ie(o,r), 

/(0,x) = /o(x), X G (0,1), 

where oj C (0,1), has been studied by P. Cannarsa, P. Martinez and J. Van- 
costenoble in [8]. Their strategy relies on appropriate Carleman estimates. 
To deal with the degeneracy at {x = 0}, they use an adequate functional 
framework that we recall in Section 2, and Hardy-type inequalities. 

The null-controllability of system (1.1) is a consequence of the internal 
null-controllability and the extension principle, since the control is located 
on {x = 1}, away from the degeneracy. The interest of the present article is 
to provide explicit controls. 

In the case of a control located on {x = 0}, an approximate controlla¬ 
bility result for a G [0,1) has been proven by P. Cannarsa, J. Tort and 
M. Yamamoto in [10] using Carleman estimates. The exact controllability 
was later proven by M. Gueye in [13] again in the weakly degenerate case 
a G [0,1) by using the transmutation method. 
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Other related one-dimensional problems have been treated; see [6, 7, 2], 
see [5] for a non-divergence setting, see [20] for a system with a singular 
potential. A multi-dimensional case has been studied in [9]. 

1.2.2. The flatness method. The main interest of the flatness method is to 
provide explicit controls. It has been developed for finite-dimensional sys¬ 
tems (see [12]) and then generalised to some infinite-dimensional systems; see 

[17] for the heat equation on a cylindrical domain with boundary control, 

[18] for one-dimensional parabolic equations with varying coefficients and 

[19] for the one-dimensional Schrodinger equation. However, the strongly 
degenerate case a G [1, 2) considered in Theorem 1.3 does not belong to the 
class concerned in [18]. Our goal is to adapt the flatness method to this 
case. 

1.3. Open questions and perspectives. The flatness method may also 
be successful on similar equations, for instance in non-divergence form as in 
[5]. For the time being, this is an open problem. 

1.4. Structure of the article. In Section 2 we recall a well-posedness 
result and the functional framework. In Section 3 we derive, thanks to 
an heuristic method, an explicit solution of system (1.1) consisting on a 
formal series development. We prove its convergence, provided that the 
corresponding flat output is in a Gevrey class. In Section 4 we discuss 
the spectral analysis of the associated stationary problem. In Section 5 we 
study the regularising effect of system (1.1) when u = 0. In Section 6 we 
construct an appropriate flat output steering the solution of (1.1) to zero, 
which concludes the proof of Theorem 1.3. Finally, we give in Appendices 
A and B a brief account of some results concerning the Gamma and Bessel 
functions needed in the proofs. 

1.5. Notation. Since all the functions appearing in the article are real¬ 
valued, we omit any explicit mention by writing, for instance, T^(0,1) in¬ 
stead of L^((0,1);M). If h G ^^([ti,t 2 ])j for some ti,t 2 £ 1^ with ti < t 2 
and A; G N*, we will denote by h'(t) and h"{t) its first and second derivatives 
and by h^'^^t), for every n G N, 2 < n < /c, the n—th derivative. 

If /ii, ^2 : M —)• R are two real-valued functions and /r G R, we will write 
/ii ~ /i 2 as X —^ to denote that lim^^^ = 1. 

We will denote by (•, •) the inner product in T^(0,1). 


2. Well-posedness 

We consider, for T > 0 and /o G T^(0,1), the following system 

dtf{t, x) - dx (x“(9a;) /(t, x) = h{t, x), {t, x) G (0, T) x (0,1), 


(2.5) 


(x“9^)/(t,x)|a;=0 = 0, 
/(Tl) = 0 , 

/(0,x) = /o(x). 


t G (0, r), 
t G (0,r), 
X G (0,1). 
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We recall below a well-posedness result for system (2.5) proven originally 
in [7]. The strategy of the proof consists in a semigroup approach and the 
introduction of adequate weighted Sobolev spaces, that we recall below. We 
refer to [7, 4] for further details. 

We introduce the weighted Sobolev space 

77^(0, !):={/£ T^(0,1); / is loc. absolutely continuous on (0,1], 

X?/' G ^2(0,1) and /(I) = o| , 

endowed with the norm 

11/11^1(0,1) •= Il/lli2(0,l) + 11 *^/^IIl2(0,1)’ ^/ 

We remark that 77^(0,1) is a Hilbert space with the scalar product 

(2-6) (/,5r)Hi := / f{x)g{x)dx+ [ x°^ f {x)g'{x) dx, V/, 5 G 77^(0,1). 
Jo Jo 

PROPOSITION 2.1 ([7], Proposition 3.2 and Theorem 3.1). Let 
/ D{A) := {/ G 771(0,1); x“/' G 77^(0,1)} , 

^ ^ 1 ^/:=-(x“/')'. 

Then, A : D{A) —>• L^(0,1) is a closed self-adjoint positive operator with 
dense domain. As a consequence, A is the infinitesimal generator of a 
strongly continuous semigroup, and for any fo G L^(0,1), andh G L^((0, T)x 
(0,1)) there exists a unique weak solution of system (2.5), i.e., a function 
f G ^0([0,r];L2(0,l)) nL2(0,T; 771 ( 0 , 1 )) such that 

f{t) = S{t)fo+ f S{t — s)h{s) ds, mL^(0,1), V7 g[0,T]. 

Jo 


As a consequence, using classical arguments (see for instance [11, Section 
2.5.3]), we deduce the following result. 

COROLLARY 2.2. Let T > O, /o G L^(0,1) and u G 77i(0,T). Then, 
system (1.1) has a unique weak solution (see Definition 1.2). 


Proof. Let /o G 7y^(0,1), u G 77i(0,T) and 

6{x) := x^, X G [0,1]. 


We consider the system 

{dt-da:{x°^daO)g{t,x) = Id{t,x), 
{x°^dfi)g{t,x)\x=o = d, 

5(7,1) = 0 , 

5 ( 0 ,x) = /o(x) - u(0)6i(x). 


(i,x) G (0,r) X (0,1), 
t G (0, T), 
t G (0, T), 

X G (0,1), 


with 


H{t,x) \= —u'(t)6{x) — u{t)A9{x), V(t,x) G (0,T) x (0,1). 
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Since H G L^((0,T) x (0,1)), by Proposition 2.1 there exists a unique weak 
solution g G ^°([0, T]; L^(0,1)) H L^(0, T; 1)) of this problem. We set 

f{t,x) := g{t,x) +u{t)e{x). 

Then, using the integral formulation associated to g, one shows that / is a 
weak solution of system (1.1) in the sense of Definition 1.2. 

The uniqueness follows since, if /i and /2 are weak solutions of (1.1), then 
/i ~ /2 is the unique weak solution of system (2.5) with h = 0, and then by 
Proposition 2.1, /i — /2 = 0. □ 


3. Explicit solution 

3.1. Heuristics. We consider the following formal expansion 

oo 2 ^ 

= '^C2k{t) , V(t,x) G (0,T) X (0,1). 

k=0 

where {c 2 k{t))keN is a sequence of real numbers. We formally have 

a — 1 


d^{x^d^f){t,x) = ^C2(fc+i)(t)(2-a)^(A: + l) 


k=0 

oo 


fc + l + 


2 — a 


1 -- 
X 2 


2k 


dtf{t,x) = ^C 2 fc(t) 
fc =0 


1 -- 
X 2 


2k 


If / solves (1.1), then the following recurrence relation holds 

C2fc(^) 


C 2 (fc+ 1)(0 — 


i2-a)^{k + l)(k + l + ^ 


yk G N. 


Choosing a flat output co(t) := y{t) and iterating, we readily have 


C2k{t) = 


(2-a)2^fe!n;=r(i + f^ 


a—1 


2-a 


, Vt G (0,T), VA; G N. 


This gives a formal solution of (1.1), 




(x^ 2 


2k 


(3.8) 

^o(2-a)2fcfc!n-=i(j + fE^ 

and a control given by u{t) = /(t, 1), which is 

y(^){t) 


u{t) = 


=0 (2 - a)2*^fc! n =1 {j + §=i 


(3.9) 
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3.2. Pointwise solutions. The goal of this section is to introduce a notion 
of pointwise solution of system (1.1) to give a sense to the heuristics made 
in the previous section. 

We define 

‘^q(0, 1) ;= {/ G ^^{[0, 1]) n ^^((0,1)) such that x°'f'{x) G ‘^°([0,1))} . 

DEFINITION 3.1 (Pointwise solution). Let ti,t 2 G M with ti < t 2 - Let 
fti G ^^(0,1) and u G ^°([ti,t 2 ])- ^ pointwise solution of system 

dtf{t,x) - dx{x‘^dxf)it,x) = 0, (t,x) G (ti,t 2 ) X (0,1), 

iqin'i x^dxfit,x)\x=o = 0, t€{ti,t2), 

f{t,l) = n{t), te{h,t 2 ), 

= /ti(x), X G (0,1), 

is a function f G ^°([ti,t 2 ] X [0,1]) n^^((ti,t 2 ) X (0,1)) such that 

(1) ■) £ ^a(O) 1): G (tl,t2): 

(2) dtf - dx{x°‘dxf) = 0 pointwisely in (^ 1 ,^ 2 ) x (0,1), 

(3) \\va.^^Q+ x^dxf{t,x) = 0, Vt G (ti,t2), 

(4) /(t,l) = u{t), yt G {ti,t2), 

(5) /(ti,x) =/ii(3:), Vx G (0,1). 

REMARK 3.2. The usual energy argument proves that, givenu G ^^{[ti,t 2 ]), 
the pointwise solution of system (3.10) is unique. We observe that, ehanging 
parameters adequately in Definition 1.2 a pointwise solution of (3.10) is also 
a weak solution. 


3.3. Convergence. The goal of this section is the proof of the following 
result. 


PROPOSITION 3.3. Let ti,t 2 G M, with ti < t 2 . If y ^ ^^([ti,t 2 ]) for 
some s G (0, 2), then 

(1) the control u given by (3.9) is well defined and belongs to i#®([ti, ^ 2 ]), 

(2) the function given by (3.8) is a pointwise solution (see Definition 
3.1) of system (3.10) in (ti,t 2 ) x (0,1) with u given by (3.9) and 
initial datum 


ft Ax) ■■= 


00 




1 -- 
X 2 


2k 


K(2-a)“(:!nL(i+e' 


Vx G [0,1]. 


Proof. Let M, i? > 0 be such that ^ n. G N, t G [ti,t 2 ]- 

Step 1 : We prove that u is well defined and belongs to ([ti, ^ 2 ]) • 
For any t G [ti,t 2 ], A: G N*, we have, as ^ 5 ^ > g, 

|?/(^)(t)| ^ Mk\^ _ M 

(2 — a)'^^k\ 1 ij + ~ R^{‘2‘ — R^{2 — 
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Hence, the series in (3.9) converges uniformly w.r.t. t e [^ 1 ,^ 2 ] and 
u G t 2 ])- Furthermore, for any n G N*, the function Cn,k{t) ■= 

- — - — - satisfies 




\Cn,k{i)\ < 


M{k + ny/ 

i?n+fc(2 _ 


Mt G [ti,t 2 ], k,n gN. 


Thus, Ylk^n,k{t) converges uniformly w.r.t t G [^ 1 ,^ 2 ]- Whence, u G 
and for every n G N, t G [ti,t 2 ], = Yl'^=oCn,kit)- 

Step 2 : We prove that u is Gevrey of order s. 

Let n G N. We deduce from last inequality that 




. 11 ) 


CX) 

s E 

fc =0 

< M 


M{k + ny.^ 

Rn+k(^2 - af'^k^ 


00 ^ 


Lfc=0 


(A:!)2-* \R{2-aY 


s\n 


R 


nl 


where we have used (A.41). The D’Alembert criterium for entire 
series shows that, whenever s G (0,2), the series above converges, 
which shows that u G ). 

Step 3 : We show that the function / given by (3.8) is well defined and 

/G^°([tl,t2]x [0,l])n^l((tl,t2)x ( 0 , 1 )). 

Let, for every A: G N, 


fk{t,x) 


(2 - a)2fcA:! Y\']=i (i + fE^) 


V(t,x) G [ti,t2] X [0,1]. 


Then, 

M / 1 \ ^ 

\fk{t,x)\ < [ R(^2-a) ) ’ e X [0,1]. 

This proves that '^^fk converges uniformly w.r.t. {t,x) G [^ 1 ,^ 2 ] x 
[0,1]. Thus,/G^0([ti,t2]x [0,1]). 

We observe that 3feo = ko{a) G N* such that {2 — a) ko > 1. Then, 
for every k > ko, /^(t, •) G ^^{[0, 1]) and 


y(^)(t)2A;(l-f)x-t 

2 ^ 

2fc-l 

(2-a)2U!n*.l( 

i + lG 



- 2) k\^-^ {r{ 2 - ay) ’ 

since (l — {2k — 1) — f > 0. This proves that Ylk>ko ^xfk con¬ 

verges uniformly w.r.t. {t,x) G [ti,t 2 ] x [0,1]. Thus, f{t,-) G 
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'^^((0,1]) for every t G [ti,t 2 ]- Note that / may not be differ¬ 
entiable w.r.t. X at X = 0 because of the hnite number of terms 
Y!k=odxfk- Moreover, dxfit,x) = for every (t,x) G 

(ti,t2) X (0,1). 

A similar argument shows that, for every x G (0,1), fi-,x) G 
^ 2 ) and 

CX) 

(3.12) dtfit, x) = 9tfk{t, x), V(t, x) G (ti, ^ 2 ) X (0,1). 

k=0 

Finally, since the partial derivatives of / exist and are continuous in 
{ti,t2) X (0,1), / G ^^{{ti,t2) X (0,1)). 

Step 4: We show that f{t, •) G ^^(0,1), for every t G (ti,t 2 )- 

Let ki = ki{a) G N* such that ki{2 — a) > 2. Working as in 
Step 3, we see that Yl,k>ki ^xfk converges uniformly w.r.t. (t,x) G 
(ti,t 2 ) X (0,1). Thus, f{t,-) G *^^(0,1), Vt G {ti,t 2 ). Furthermore, 

2(fc-l) 




1 -- 
X 2 


(3.13) d^{x'^d^f){t,x) = J2 - , , 

(2 - a)2(fc-b(fc - 1)! ^ 

for every {t,x) G (ti,t 2 ) x (0,1). From Step 3, we obtain 

y^^\t) 2 k (1 - f) 


\x^dxf{t,x)\ = 


y 

S {2-aY^k\\{y[j + ^ 


< 2M 1 - 


a 


E 

k=l L 


k 


k'?-^ \R{2-aY 


X, 


for all (t,x) G (^ 1 ,^ 2 ) x (0,1), which implies, since a G [1,2), that 

x^'dxfit, x) —)• 0, as X —O’*". 

Therefore, f{t, ■) G for every t G (ti,t 2 )- 
Step 5: According to (3.12) and (3.13), an straightforward computa¬ 
tion shows that the equation in (3.10) is satisfied. 

□ 

4. Spectral Analysis 

The goal of this section is to give the explicit expression of the eigenfunc¬ 
tions and eigenvalues of the spectral problem 

Aip{x) = \ip{x), xG(0 , 1), 

{x°^^') |x=0 = <^(1) = 0, 

where A is given by (2.7). We will make use of several results about Bessel 
functions recalled in Appendix B. Form now on, we use the notation 

a — 1 
2 — a 


(4.14) 


(4.15) 


v \= 
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PROPOSITION 4.1. Let 

(4.16) ipk{x) = I S (0> 1)> k gN*. 

Then, 

(1) G D{A), Vfc e 

(2) (pk satisfies ( 4 - 14 ) with 

(4.17) Afc:=(l-|)'4„ Vfeer, 

(3) {pk)keN* ® Hilbert basis o/L^(0,1), 

(4) for every Jq G L^(0, 1) the solution of (2.5) with h = 0 writes 

OO 

(4.18) f{t) = Y,e-^’‘\fo,Tk)Tk inL\0,l),ytG[0,T]. 

k=l 


Proof. We will note for simplicity bk ■= and (pk ■= ^Tk^ for every 

kGW. 


Step 1: We prove that pk € D{A), for every kGW and that Apk — 
^kPk — 0 . 

Let kGW. We observe that pk G ‘^°°((0,1]) fl ^°([0,1]), for any 
kGW and x G (0,1). We have 


(4.19) p'kix) 


1 — 1 + a , 1 — “ \ 

^ X 2 Jjp^Ji/^kX Ju^k 



°‘J({ju,kX 



Whence, using (B.48) and Lemma B.3, we deduce 




(1 — a) O (x^ + O 

X— >-0+ V / x^0+ 



(4.20) 


It follows that x'ip'n G L^(0,1). Thus pk G R^(0,1). Moreover, 
from (4.19), a direct computation shows 



Then, evaluating equation (B.47) at z = ju^kX^ 2 and multiplying 

Q: — 3 

by x~^, it follows 
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Substituting in (4.20), this gives 

Then, we readily have € /4^(0,1) C L^(0,1). Thus, ifk S 

D{A). Moreover, Aipk = Xk^k- 

Step 2 : We check the boundary condition of (4.14) at x = 0. 

We observe first that the case a = 1 is straightforward. From 
(4.19), (B.48) and Lemma B.3, we have 

\x'"(p'n{x)\= O . 

Then, it follows that \\m.^^Q+ x'^(p'^{x) = 0. This shows, combined 
with Step 1, that satisfies (4.14). 

Step 3 : We prove that {'Pk)k£n* is an orthonormal family in L^(0, 1). 
Let n, m G N*. Then, changing variables and using (B.46), we get 


/ ^n{x)Pm{x)dx 

Jo 

l—a ^i-'{ji-',nX 2 ) Jyi^jn^YnX 2 ) 


= ( 2 -«: 


/ 


dx 


\Jv+l{jy,n)\ \Jv+l{jii,m)\ 

fl 

yJuijujny^Juijuymy) dy = 6n 


\Ju+l{jiy,n)\\Ju+l{jiy,m)\ JO 
where 5n,m stands for the Kronecker delta. 

Step 4: We prove that {pk)ke'^* is a Hilbert basis of L^(0, 1) by check¬ 
ing the Bessel equality. Let / G T^(0, 1) and let 


(4.21) 


CLk • — 


:= / f{x)pk{- 

Jo 


x)dx, VfcGN*. 


Then, using Lemma B.l and changing variables twice, we get 


oo oo „ 

k=l k=l 


\/2 — a i-a 


/O \Ju+l[ju,k)\ 


X 2 Ju [ju,kX^ 2 j dx 


2 — a 


E 

k=l 


-1 


V2y 


/ 2/2-« 2/(y2-a) 

'0 |4y+UJjy,fcj| 


Jv{ju,ky) dy 


2 — a 


/■l 2(n-l) 

/ y 

Jo 


+1 


fiy 


2 

2-a ' 


dy 


\f{z)fdz=\\f\\l2^o^^y 


Step 5: Finally, (4.18) is a consequence of [3, Theorem 8.2.3, pp.237” 
240]. 

□ 
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5. Regularising effect 


We use the orthonormal basis obtained in Proposition 4.1 and some prop¬ 
erties of Bessel functions to quantify the smoothing of the solution of system 
(1.1) when u = 0. 

PROPOSITIONS.!. Letfo G 1.2(0,1), T > 0 and let f £ '^“([0, T]; L2(0, 1)) 
he the unique weak solution of system (2.5) when h = 0, aceording to Propo¬ 
sition 2.1. Then, there exists Y G ‘^°°((0,T]) such that for every a G (0,T), 




and 


(5.22) fit, x) = Y, ^ ^ 


2n 


^o(2-a)2-n!m=i(j + fE' 


V(t,x) G [a,T] X [0,1]. 


Moreover, f solves system (3.10) pointwisely (see Definition 3.1) in (ct, T) x 
(0,1) with u = 0 and initial datum /o-(x) = /(a,x). 

Proof. Let u be given by (4.15) and as in (4.21). Let a G (0,T) be fixed 
but arbitrary. Let t G [cr,T] be fixed. By (4.18) and (B.43), we have, for 
a.e. X G [0,1], 


/(M) = .yj. 


k=l 

oo 


u I 3v,nX 2 




, e 
k=l 

OO OO 


i-iy 


' . l-£L 

Jl2,kX 2 


2n-\-u 


(5.23) = J2J2Bn,kit,x), 

k=l n=0 

where, for every (n, /c) G N x N*, 


nir (n + l + i/) 

n=0 


n -271+1/ 


(—lYi 

Bn,k{t,x) ■.= e~^^%k 


1 -- 
X 2 


2n 


n!r(n + l + u)22^+^^|J^+i(j^_fc)|’ 


and bk := ak^/2 — a, \/k G N* 
Step 1 : We show that 


CXD / OO 


(5.24) 


X] ( X] Va;G[0,1]. 

n=0 Vfc=l / 

Indeed, since > 0, we have for every (n, A:) E N x N* and x E [0,1], 


I ^n,k (Aj X^ I ^ 


\h\j, 


2n+u -Xkcr 
u,k ^ 


22^+^n!r(n + l + u)|J^+i(j^,fc)| 
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(5.25) 


< 


^i\Ok\e 

22"-n!r(n + 1 + v) ’ 


for a constant Ci > 0, using Lemma B.4. 

We fix n G N and we define the function /i“ G M+) by 

/i"(x) := , Vx G [0, +oo), 

which satisfies that 

(5.26) -^K{x) > 0, Vx G (0, iV“) and -^K{x) < 0, Vx G (K,oo), 

UvO Clt/L 


where A^" := 2 ^J a 4 ( 2 -a) )' Hence, from (5.25) and (4.17), 


(5.27) 


^ \Bn,k{t,x)\ < 


Cl supfc \bk\ 


'^K{ju,k) 


2‘^^n\r(n + 1 + i^) ^ 

k=i ^ ' k=l 

Introducing iL“ := sup{/c G N*]ju,k < we write 


(5.28) K{jy,k) = KUu,k^) + hnUu,K^+i) + Y1 K{jv,k) 

fceN*-{A'“,A'“+l} 


k=l 

On one hand, we have 


K{3 u,k^) + K{j^,K^+i) < ‘^KiK) 


< 2e r+4(2-0) 7 I 


a 


4(2 - a) 


n+ 


4(2-0) 


cr V 2 — a 


-I n+ 


4(2-o) 


(5.29) < C 2 r n + 


a 


1 

4(2 -a)^2 


1 


cr \ 2 — a 


-I 72+ 


4(2-o) 


for a constant C 2 > 0, using Lemma A.l with a = 1, 6 = ^. On the 
other hand, using (5.26), we write 

K{ju,k) < 

fc6N*-{K“,A'“+l} 

1 /-Lx+i “ 1 /-Lx 

< Y -/ h^{x)dx+ Y -■—- / h^{x)dx 

3v,k+l Jv,k J j^a_^_i3iy,k Ju,k-l J j„ 


1 




< sup 
fcSN* L 3u,k+l Jv,k J 

roo 

< C'a / h^{x)dx, 

Jo 


K 


' 3v,k — l 

'(x)dx+ f /j“(x)dx| 
27“ 


Ji^,Kn + l 
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for a constant 6*3 > 0, using (B.45). Moreover, we have 


f 


/i“(x)dx = 




roo . 

r--i^s 

1 2n-\- 


dx 


2nH- 


2(2-a) 


1 


^/a \2 — a 


2(2-0) 


2y/^ \2 — a 
+ 1 /*oo 




/ 


4(2-0) 2 


1 


2 [\2 — aJ 


2 "-+ 2 ( 2 - 0 )+^ 


r n + 


a 1 

4(2 - a) 2 


where we have used (A.38) withp = + ^. Hence, combining 

this with (5.28) and (5.29), we get 

2n+ 


J2K{ju,k)<(C2 + 


k=l 


C 3 


v ^(2 - a) 


— (— 

y/a \2 — a 


■ 2 ( 2 - 0 ) 


r n + 


a 


1 

4(2 -a)^2 


which, according to (5.27), implies 


^ \Bn,k{i,x)\ < C4 


k=l 


[ 1 



V 2 - a)_ 


2 ’^+ 2 ( 2 ^ T (n+ 4 ( 2 “^) + 

22 ’^nir (n + u + 1 ) 


Henceforth, the D’Alembert criterium for entire series gives (5.24). 

Step 2: We find Y e '^^{[a,T]) such that (5.22) holds. 

Thanks to Fubini’s theorem, (5.23) and (A.39), we may write 

“ A(2_c)2»„!n".i0 + o’ 


yn(t) := 


n=0 

where, for every n € N, 

(l)n^ 2 ^ (1 a) 


2 n (X) '272+1^ 


2 ^r 


1 

2-a 


E 

/C=l 


dp^c 


\Ju+l{jii,k)\ 


Vte [a,T], 


and u is given by (4.15). Putting 


(5,30) Yit) := t e ^r], 

2"r (ji;) St l■'■'+lG.‘)l 

we have that, since u > 0, T is analytic in [u, T]. Moreover, 

Y^^\t) = yn{t), VteKr],VnGN. 

Hence, we obtain (5.22) with this choice. Since cr G (0, T) is arbi¬ 
trary, we have in addition that Y G ^°°((0,T]). 
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Furthermore, applying Proposition 3.3 to (5.22) with ti = a and 
t 2 = r, we deduce that / solves (1.1) pointwisely in (cr, T) x (0,1) 
with M = 0 and fcr{x) = f{a,x). 

□ 


6. Construction of the control 


Let s G M with s > 1. The function (see [17, Section 2] and [21, Theorem 
11.2, p.48]) 


(6.31) 


Mt) 


1 , 

-- 1 ) 

e-(i-t) ^+e-t ^ 

0 , 


if t < 0, 
if 0 < t < 1, 
if t > 1, 


belongs to 1^®([0,1]) and satisfies 

(6.32) ,^,(0) = 1, <^,(1) = 0, #(0) = #(l)=0,ViGNT 


Proof of Theorem 1.3. Let /o G L^(0,1), T > 0. Let / and Y be given by 
Proposition 5.1. 

We pick r G (0, T), s G (1, 2) and we set the flat output 


y{t) ■= 4>s 



VtG (o,r]. 


which belongs to ^°°(0,T). Moreover, for every a G (0,T), y G ?f^([(T,T]), 
as it is a product of two functions in l^®([cj, T]). We define accordingly the 
function 




E 




2n 


T^i(2-a)2-n!n"=i(j + f^' 


2-a 


V(Lx) G (0,T] X [0,1], 


and the control 
(6.33) 


u{t) 


0, iG[0,r], 

/(t,l), fG(r,r]. 


Since y G ?f^([cr, Tj) for some s G (1,2), Proposition 3.3 shows that 


(6.34) 


Vcr G (0,T), / is the pointwise solution of (3.10) with 
ti = a,t 2 = T, /tj = f{a, •) and (6.33). 


As a consequence of (6.32), we have 


y(t) = y(t), Vt G (0,r], 

(6.35) y{T) = 0. 

Whence, f{t,x) = f{t,x), for every {t,x) G (0,r) x (0,1). Thus, as / G 
‘^‘^([0, Tj; L^(0,1)), we deduce 

(6.36) /G'^°([0,r];L2(o,l)), 

(6.37) /(0)=/o inL2(0,l). 
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We have to check that / is the weak solution of system (1.1) on (0,T). To 
do so, and according to Definition 1.2, let t' € (0,T) and let -0 satisfying 
(1.2) and (1.3). Then, by (6.34) and since a pointwise solution is a weak 
solution (see Remark 3.2), we have, for every cr > 0, 


[ [ f{t,x){dt'ilj + da;{x°‘dxip)){t,x)dtdx 

Ja Jo 

— [ f{t',x)'ip{t',x)dx— f f{a,x)'ip{a,x)dx+ f u{t) (x^dxip) {t,l) dt. 
Jo Jo Jo- 

Then, from (6.33), (6.36), (6.37) and (1.2), taking cr ^ O’*", we get the 
conclusion. 

Finally, by construction (6.35) implies that f{T, x) = 0, for every x G 

( 0 , 1 ). 

□ 


Acknowledgements. I thank Karine Beauchard for suggesting me this 
problem and for many fruitful discussions. 


Appendix A. Some properties of the Gamma function 
For any p G M"*", the Gamma function is defined (see [1, 6.1.1, p.254]) by 

rOQ 

(A. 38 ) r(p) := / dt, 

Jo 

which is a monotone increasing function on (0, oo). Furthermore, (see [1, 
6.1.15, p.256]) 

(A.39) F(x + 1) = 3:r(3:), Va:G(0,oo). 

We have the following asymptotics of the Gamma function. 

LEMMA A.l ([l], 6.1.39 ). Let a G and 6 G M. Then, 

(A.40) T(ax + b) ~ 

ai—)-oo 


We show an inequality used in Proposition 3.3. 

LEMMA A.2. 

(A.41) {n + k)\ <2'^+^n\k\, Vn, A; G N. 

Proof. Let us observe first that 

(A.42) (2n)! < 2‘^^nl^, Vn G N. 

This inequality follows by induction, since, for every n G N, 
(2(n + l))! = (2n)!(2n +l)(2n + 2) 

< (2n)!2^(n +1)2 < 2^(’"+^)(n +1)!. 
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To show (A.41), we assume, w.l.o.g., that n < k. Then, using (A.42), 

k—n 

{n + k)\ = (2re)! JJ(2n + j) 

1=1 

k—n 

< {2n)l2^-^ n 
1=1 

□ 


Appendix B. Some properties of Bessel functions 

Let p £ M. The Bessel function of order u of the first kind is ([1, 9.1.10, 
p.360]) 

{B.43) JJ.) := (£)“*-, V. € [O.oo). 

n!l n + u + 1 V2/ 

n=0 ^ ^ 

We denote by the increasing sequence of zeros of J^, which are 

real for any u > 0 and enjoy the following properties (see [1, 9.5.2, p.370] 
and [15, Proposition 7.8, p.l35]). 

(B.44) U < jii^n < Ji/,n+l) G hi , 

(B.45) ju,n+i - ju,n -5- vr, as n oo. 

We also have the integral formula ([1, 11.4.5, p.485]) 

1 

(B.46) / yJu{ju,ny)Ju{ju,niy)dy = -\Ji,+i{j,y^n)f6n,m, Vn,mGN*. 

Jo 2. 

This allows to show the following. 

LEMMA B.l. [14, p.40] Let u > 0. The family {ri’nlngn* defined by 

\f2z 

Wn{z) ■= T-; 'dz G (0, 1), 

I <-' 1^+1 \Ju,n)\ 

is an orthonormal basis of L^{0, 1). In particular, if f G L^(0,1) and dn := 
fo fiz)wniz) dz, Vn G N*, then ||/||i 2 (o,i) = Er=i l^nP- 

We recall that Vu G M, the Bessel function satishes the following 
differential equation (see [1, 9.1.1, p.358]) 

(B.47) z^J"{z) + zJl{z) + (z^ - u^) Jy(z) = 0, Vz G (0, Too), 

and the recurrence relation (see [1, 9.1.27, p.361]), 

(B.48) 2J;[,(2;) = Jv-i{z) + Jy+i{z), Vz G (0, Too). 
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Asymptotic behaviour. We recall the asymptotic behaviour of for 
large arguments and near zero. 

LEMMA B.2. [15, Lemma 7.2, p.l29] For any u G M, 

’ y irz \ 2 i) i-tooGViy 


LEMMA B.3. [1, 9.1.7, p.360] For any u G M \ {— N*}, 


Ju{z) 


z^Q 2^r(u + !)■ 


The following asymptotic result is important in the proof of Proposition 
5.1. We give the proof for the sake of completeness. 

LEMMA B.4. Let v G M"''. Then, 


(B.49) \/ jv,k\Ju+l{ju,k)\ — \ ^ ^ 1 ■ 

V vr k^co \Ju,k J 

In particular, there exists a constant Ci > 0 such that for all k gW , 


1 

Jy+l {ju,k)\ 


— Cl \/ ju,k- 


Proof. Using Lemma B.2, for z/ + 1 and x = j^^k, 


COS 

sin (jjyh - 


2 

7TU TT 


+ O , . 

fc^-OO \Ji^k 


1 


+ o 

k^OO \jyk 


2 4 

Using again Lemma B.2 with v and x = jy^k-, we have that 


cos ^ = O (—), 

2 4/ fc^OO \Jy^k, 


which gives 


vru TT 


Sm [Jy^k 


and then (B.49). 


\ 


1 + O 

k^OQ 


, iv,k , 


= 1+ o ( — 

k^OO \Jy^k 


□ 
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